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A light impurity in a Fermi sea undergoes a transition from a polaron to a molecule for increasing 
interaction. We develop a new method to compute the spectral functions of the polaron and molecule 
in a unified framework based on the functional renormalization group with full self-energy feedback. 
We discuss the energy spectra and decay widths of the attractive and repulsive polaron branches 
as well as the molecular bound state, and confirm the scaling of the excited state decay rate near 
the transition. The quasi-particle weight of the polaron shifts from the attractive to the repulsive 
branch across the transition, while the molecular bound state has a very small residue characteristic 
for a composite particle. We propose an experimental procedure to measure the repulsive branch in 
a 6 Li Fermi gas using rf-spectroscopy and calculate the corresponding spectra. 
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I. INTRODUCTION 

A single impurity 4,-atom immersed in a background of 
f-fermions is a screened Fermi polaron 1J below a critical 
interaction strength. For larger interaction the ground 
state changes its character and the impurity forms a 
molecular bound state with the bath atoms. This qual- 
itative change marks the polaron-to-molecule transition 
and was predicted by Prokof'ev and Svistunov |[2 . The 
transition can be regarded as the modification of the 
two-body problem due to medium effects: For two non- 
relativistic particles with an interaction characterized by 
the s-wave scattering length a, a weakly bound molecule 
is the ground state for positive a, whereas for negative a 
the molecular state ceases to exist and the ground state 
is given by the free atoms. 

Related impurity models have been studied for many 
years, in particular the Kondo effect for a fixed impurity 
with discrete energy levels immersed in a Fermi sea of 
conduction electrons 3 . A mobile but very heavy impu- 
rity loses its quasi-particle character in low dimensions 
d = 1, related to the orthogonality catastrophe 4 . The 
transition of the light mobile impurity poses a new and 
challenging many-body problem. 

The polaron-to-molecule transition is observable with 
ultracold atoms where the scattering length can be tuned 
via Feshbach resonances |5j. Using radio-frequency (rf) 
spectroscopy the lineshape, ground-state energy, and po- 
laron quasi-particle weight have been measured across 
the transition [5J. The ground-state properties near the 
transition have been calculated using variational wave- 
functions [7Ul3j. non-self-consistent T-matrix approxi- 
mations [8], l/N expansions [14], Wilsonian renormal- 
ization group [15], variational Monte Carlo [TB] and dia- 
grammatic Monte Carlo (diagMC) [2]. 

It turns out that even beyond the critical interaction 
strength the polaron remains as a long-lived excitation 
above the molecular ground state, and conversely the 
molecule becomes an excited state on the polaronic side 
of the transition. The proper description of these ex- 



cited states and in particular their finite lifetime remains 
a difficult problem. Using a phenomenological model in a 
three-loop calculation, Bruun and Massignan have shown 
that the decay rates of the excited states decrease rapidly 
as Aw 9 / 2 toward the transition, where Aw is the energy 
difference between the ground and excited states |17) . 

Recently, interest has focused on an additional fea- 
ture present in the polaron-to-molecule transition: The 
renormalization of the ^.-spectral function by the strong 
interactions leads to the appearance of an additional 
quasi-particle excitation for positive energies. This ex- 
citation corresponds to a Fermi polaron interacting re- 
pulsively with the f-Fermi sea. The repulsive polaron 
has a finite lifetime which becomes very small at uni- 
tarity, a — > oo. In the opposite limit of weak coupling 
a — > 0, the repulsive polaron is long-lived and its spec- 
tral weight approaches unity. This justifies perturba- 
tive methods which neglect the presence of the molecular 
channel jT5J Qj|] . The repulsive polaron has been stud- 
ied theoretically in the context of ultracold gases by Cui 
and Zhai |20j and more recently by Massignan and Bruun 
[21j . Experimentally, the repulsive polaron has not yet 
been observed in the strong-coupling regime. However, 
it has important implications for the stability of a fer- 
romagnetic phase in Fermi gases (3U] [22H2BJ. Further- 
more, a single 4,-fermion in an f-Fermi sea corresponds 
to a strongly imbalanced, two-component Fermi gas near 
full polarization, and the polaron-to-molecule transition 
sheds light on a region of the zero-temperature phase di- 
agram of the polarized Fermi gas [101 El l2Z] • 

So far, there is a lack of theoretical work which de- 
scribes all of these features within one unified approach. 
In this paper we present a method capable of doing this: 
Using a novel numerical implementation of the functional 
renormalization group (fRG) we are able to determine 
both the full spectral functions and the quasi-particle 
features of the polaron-to-molecule transition. 

With the renormalization group one usually limits one- 
self to a few running couplings. In this work we develop 
a computational tool which allows us to keep track of 
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the renormalization group flow of fully momentum and 
frequency dependent Matsubara Green's functions based 
on an exact renormalization group equation |28) . By in- 
troducing an auxiliary bosonic field which mediates the 
interaction between f- and 4,-atoms, we are able to accu- 
rately capture the momentum and frequency dependence 
of both vertex functions and propagators. In particular 
for non-universal quantities our approach, which com- 
plements the proposal by Blaizot, Mendez-Galain, and 
Wschebor [29], may prove vital, and we demonstrate its 
efficiency in the polaron problem. 

Our main results from the fRG are the following. We 
give the spectral functions of the 4,-fermion and — for the 
first time — also of the molecule in a wide range of in- 
teraction values. From this we extract the quasi-particle 
properties including energies, decay rates and residues 
for both the polaron and the molecule. We connect the 
quasi-particle weight of the molecule with a measure of its 
compositeness introduced by Weinberg [3U]. We propose 
a new experimental procedure to measure the repulsive 
polaron in a 6 Li Fermi gas using rf-spectroscopy, and we 
employ our spectral functions to predict the expected rf- 
response. 

The paper is organized as follows: in section [IT] we in- 
troduce the model of the polaron-to-molecule transition 
and derive appropriate fRG flow equations. These equa- 
tions are solved in section [TTT] for the full frequency and 
momentum dependent spectral functions, which are pre- 
sented and interpreted in section IV Section[V]is devoted 
to rf-spectroscopy, and we conclude with a discussion in 
section [Vl] In the appendix the fRG flow equations are 
solved using the simpler derivative expansion, which al- 
ready gives qualitatively correct results. 



II. MODEL AND RG FLOW EQUATIONS 

In this work we study a two-component Fermi gas in 
the limit of extreme population imbalance at T = 0. The 
microscopic action describing the system is 

s = I { E ^ - A - + w-M^vt} (i) 

in natural units h = 2m = 1 and with imaginary time 
t. The Grassmann-valued, fermionic fields ip* and ip\ 
denote the up and down fermions, which have equal 
mass m. The associated chemical potentials \i a are ad- 
justed such that the f-fermions have a finite density 
= kp/(6ir 2 ) while there is only a single impurity J.- 
fermion. The atoms interact via a contact interaction 
with coupling constant g which is related to the s-wave 
scattering length a via g = 8ira. The T-matrix acquires a 
complicated frequency and momentum dependence in the 
strong-coupling limit. It is then convenient to perform a 
Hubbard-Stratonovich transformation of the action ([!]) 
introducing a bosonic molecule (pairing) field <j) ~ i^i^Pf 
which mediates the two-particle interaction gip^ip^ip^ip^ . 



The resulting action is given by 
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+ h(^l<l) + h.c.)\ (2) 



with a real Yukawa coupling h for the conversion of two 
fermions into a molecule. Integrating out the bosonic 
field 4> shows that ^ is equivalent to the single-channel 
model ([I]) provided that —h^G^^ = 9 an d h — > oo 
P31 EH] • As h 2 ~ AB this limit corresponds to a broad 
Feshbach resonance [5]. 

The physical properties can be accessed via Green's 
functions which are derivable from generating function- 
als. The one-particle irreducible vertex functions 
are obtained from the effective quantum action T, which 
can, for instance, be computed perturbatively in a loop 
expansion. As we are interested in the intrinsically non- 
perturbative regime of fermions close to a Feshbach res- 
onance where the scattering length a diverges we employ 
a different approach. T includes quantum fluctuations 
on all momentum and energy scales. The main idea of 
the functional renormalization group (fRG) is to intro- 
duce an interpolating effective flowing action which 
includes only fluctuations on momentum scales q > k 
larger than the renormalization group scale k. At the 
UV scale k = A — > oo the effective flowing action reduces 
to the microscopic action S which does not include any 
quantum corrections. In the infrared limit k — ¥ 0, 
equals the full quantum action T. 

The evolution, or flow, of Tfe with the RG scale k is 
given by the exact renormalization group equation [2 8) 



d k T k = l - STr 



f 



,(2) 



+ R k 



dkRk 



(3) 



The supertrace symbol STr denotes a loop integration 
over frequency and momentum as well as the summation 
over all fields and internal degrees of freedom, with a mi- 

(2) 

nus sign for fermions. T k is the full, field dependent 
inverse two-point Green's function at scale k, and R k is 
a regulator taking care of the successive inclusion of mo- 
mentum scales. The regulator acts by imposing a large 
mass term on the modes with momenta q < k lower than 
the cutoff scale k. R k can be chosen freely as long as 
Rk^o —> and R k -+cc —> oo. For further details on the 
fRG we refer to the literature [32HI2I , and to its applica- 
tion to the BEC-BCS crossover [35] . 

r^. is in general a functional of the fields and contains 
all possible operators of the fields allowed by the sym- 
metries. For this reason its exact calculation is usually 
impossible and one has to rely on approximations for T^. 
In this work we will use the truncation 
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FIG. 1: Diagrammatic representation of the fRG flow equa- 
tions Q for the impurity and molecule propagators. 



with Matsubara frequency w. For the momentum and fre- 
quency dependence of the fc-dependent, or flowing, prop- 
agators of the 4,-fermion Gj.^ and the boson G$. k we will 
present a simple approximation in terms of a gradient 
expansion in the appendix. In the following section |III| 
we develop a new numerical method to solve the renor- 
malization group flow of the propagators as completely 
general functions of u> and p. This will enable us to cap- 
ture decay rates and dynamic effects which is not possible 
in a simple gradient expansion. 

Within the truncation ( 4| the Yukawa coupling h is not 
renormalized which can be seen by a simple argument: 
In Eq. @ we neglect a term ijj*ij;*tp^tp^ which would be 
regenerated during the flow by particle-hole fluctuations 
[36j . Similarly, a term if)!^* 'ip-f4> for the atom-dimer in- 
teraction is neglected. Both terms would lead to a renor- 
malization of h. Due to their omission, however, there is 
no diagram generating a flow of h which is why d k h = 
and h remains independent of frequency and momentum. 
Furthermore, we neglect terms (</>*</>)™- 2 which would 
give higher-order corrections to the bosonic self-energy. 
The majority f-atoms are renormalized only by the sin- 
gle impurity ,J,-atom to order 1/iV-j-, hence one can neglect 
the renormalization of the f-atoms in the thermodynamic 
limit, and the chemical potential = ep = k F /(2m) is 
that of a free Fermi gas (we work in units where the Fermi 
momentum kp = 1) Ej- 

As Eq. ^ is a (functional) differential equation, it has 
to be supplemented with appropriate initial conditions at 
the UV scale A, which are obtained from few-body (vac- 
uum) physics. The s-wave vacuum scattering amplitude 
for the interaction of an f- and 4,-fermion with momenta 
q, — q in the center-of-mass frame is given by (q = |q|) 
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— 1/a — iq 



(5) 



f(q) is related to the full molecule propagator G™r eval- 
uated at the infrared RG scale k = 0, 
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GS(^2g 2 ,p=0), 



(6) 



where u> = 2q 2 is the total kinetic energy of the interact- 
ing atoms. The subscript R indicates that the analytical 
continuation to the retarded function of real frequencies 
(iu) — > u) + iO) has been performed. In [37] the exact vac- 
uum molecule propagator G™j^ has been calculated using 
the fRG and agrees with the well-known result 
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This expression for G™^ reproduces the correct scatter- 
ing amplitude Q when inserted into Eq. (|6| , and dictates 
the form of the UV propagator G^\ for a given choice of 
regulator. 

Furthermore, the initial condition for the fermions 
is given by their form in the microscopic action pj, 
G^^, =A (w, p) = — iu + p 2 fj, a . As we will see in the 
following, the momentum and frequency dependence of 
both propagators, G,p ik and Gj.^, is strongly renormal- 
ized during the fRG flow toward the infrared, which leads 
to a rich structure of the spectral functions. 

After inserting the truncation Q into the flow equa- 
tion ([3]), the flow of G^. and G^fc is derived by taking the 
appropriate functional derivatives of Eq. ([3]) with respect 
to the fields. One obtains the fRG flow equations 



dkP^uiP) = h 2 d k / G% M {Q)Gl k (Q - P) 

JQ 

d k P^, k {P) = ~h 2 d k [ Gl k (Q)Gl k (P - Q) (8) 



Q 



with the multi- indices P — (w,p) and Q = (y, q). The 
P k = GjT 1 on the left-hand side are the flowing inverse 
propagators without the regulator from Eq. Q , while the 
propagators G c k on the right-hand side are regulated: 



G k = 1/P fe 



G c k = l/(P k + R k ). 



(9) 



The tilde on d k indicates that the derivative with respect 
to the RG scale k acts only on the regulator term R k in 
the cutoff propagators G k : specifically, the single-scale 
propagators read d k G k = ~{G k ) 2 d k R k in Eq. (|8]). Note 
that the flow equations ^ , which are depicted in Fig. [I] 
have a one-loop structure but contain the full propaga- 
tors at scale k: by integrating the flow, diagrams of arbi- 
trarily high loop order are generated and constantly fed 
back into each other. It is especially for the latter rea- 
son that our approach goes beyond other approximations 
used for the description of the polaron problem such as 
for example the non-self-consistent T-matrix approxima- 
tion [8] [21] . The goal of this paper is to solve the system 
of differential flow equations Q. 

In the following we choose sharp cutoff functions R k 
which strictly cut off all momentum modes with |p| < k 
while the frequencies are not restricted |54j . Then the 
regulated Green's functions G k take the particularly sim- 
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pie form (for unoccupied 4,-atoms) 

6(]p\-k) 



G|, fc (w,p) 
P) 



Pj. )fc (w,p)' 
g(|p|-fc) 
P^,fc(w,p) ' 
fl(|p 2 - Mt |-fc 2 ) 
Pt,feO*>,p) 



(10) 



For the f-atoms it is crucial to regularize the low-energy 
modes around the Fermi energy fj^. Using the Dyson 
equation 

Pk(u,p) = G^(lj, P ) = Go X (w,p) - E fc (w,p), 

where Go = G k =A denotes the free (UV) and G k the full 
Green's function at scale k, we can in particular identify 
the gap term 



:=n, fe (0,O) = -^-£ ; , fe (0,O). 



(11) 



The chemical potential /i^ is the energy required to add 
one J,- atom to the system, 



H = P(iV+) - E(N± - 1), 



(12) 



and is independent of the cutoff scale k. The interac- 
tion effects on the ^-fermion, which are successively in- 
cluded during the flow, are captured by the fc-dependent 
self-energy Efc(w,p). In the polaron problem we are in- 
terested in a two-component Fermi gas in the limit of 
extreme population imbalance where one considers only 
a single |-atom, Ni = 1, and relation (12) is used to 



determine the ground-state energy of the system. This 
value of the chemical potential /i^ then marks the phase 
transition from a degenerate, fully polarized f-Fermi gas 
to a phase of finite j-fermion density [TU]. Accordingly, 
for all choices of /x^ < //|, there has to be a vanishing oc- 
cupation of both 4,-fermions and molecules at every RG 
scale k, which leads to the condition 

p M (o,o,/4)>o, p M (o,o lA 4)>o V/4< w . (13) 

In order to have only a single 4,-atom or molecule, 
has to be determined self-consistently such either the 
J,-atom or the molecule <f> is gapless in the infrared, 
Pi/0,fe=o(0, 0) = (ground state). We find that depend- 
ing on the value of the dimensionless coupling (fcpa) -1 , 
either the polaron or the molecule becomes the ground 
state, see Table [IJ The polaron-to-molecule transition oc- 
curs at the critical interaction strength (fci?a c ) _1 at which 

m l = m l = °- 



coupling 


(k F a) 1 < (k F a c ) 1 


(k F a) 1 > (k F a c ) 1 


ground state 
1 gap 
<t> gap 


polaron 
n(0,O) = mJ = 
P*(0,0) = m3 > 


molecule 
P±(0,0) = ml>0 
P+(Q,0) = m% = Q 



TABLE I: Conditions for the polaron and molecule ground 
states. 



III. RG FOR FULL SPECTRAL FUNCTIONS 

The main goal of this work is to solve the system of flow 
equations ^ without imposing any constraints on the 
frequency and momentum dependence of the polaron and 
molecule propagators. This problem can only be solved 
numerically, and the inverse, flowing Green's functions 
Pl/cf, ,fc(w, p) are evaluated on a discrete grid in frequency 
and momentum space, 



Pl,fc(w,p) 
Pm(w,p) 



i,k 
pi] 



(14) 



We choose a logarithmically spaced, finite grid with 
u>i £ (0, . . . , w max ) and p 3 € (0, . . . ,p max ). As a result 
of rotational invariance the propagators depend only on 
the magnitude of spatial momentum p = |p|, and due to 
the condition P*(u>) = P(— w) for Euclidean (Matsubara) 
propagators it is sufficient to consider positive frequen- 
cies only [35]. The full (oj,p) dependence is reconstructed 
from the finite number of couplings Pu^ k by cubic spline 
interpolation, 

P l/4>A U 'P) = Spline({P 4 ' 



3 

= c 1/m( w - w <) { (p- 



Pj)\ (15) 



with uj € (u>i, Ui+i), p £ (pj,Pj+i) and k the corre- 
sponding spline coefficients. For the asymptotics of the 
propagators for high frequency and momentum modes, 
ui > w max and/or p > p max , we choose the simple fit 
models 
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M4 



(16) 
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with ffak determined by a continuity condition from the 
numerical value of P^^ for the largest momenta |p| = 

Pmax- 

In order to keep the numerical cost of computing the 
flow equations Q low it is advantageous to employ the 
sharp momentum regulator functions R^k, Pf.fcj and 
P^fc defined by Eq. (p~0|) . This reduces the number of 



loop integrations by one. The flow equations evaluated 
by our algorithm are then given by 



9fePj.,fe(a;,p) 



d k P^ k (u}, p) 



h 2 
(2tt) 3 



dx / dv 

1 J-ao Jo 

Xk(p,q,x) 



q 2 dq 



P<i>.k(v, q)Pf,k(v - w, q - p) 
h 2 i- 1 



(27T) 



dx 



dv I q dq 
lo 



PiA^ q)Pr,fc( a; - ^ p - q) 



(17) 
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where we have defined the characteristic function 



IV. FULL SPECTRAL FUNCTIONS 



Xk(p,q,x) = 5(q - k)9(\{p - q) 2 - 
+ 2k6(q-k)S(\(p-q: 



Mtl 

2 



Mtl 



(18) 



and x = cos 9 expresses the angle 9 between the momen- 
tum vectors p and q such that |p±q| 2 = p 2 +q 2 ± 2pqx. 

The initial condition at the UV scale k — A is deter- 
mined by the few-body calculation ([7|. As we employ 
no approximation for the momentum and frequency de- 
pendence of the molecule propagator we are able to in- 
corporate the exact two-body scattering amplitude ([5]), 
in contrast to the calculation in the appendix using the 
derivative expansion where this is not possible. The vac- 
uum problem can be solved exactly using the sharp reg- 
ulators (10 1, which leads to the UV molecule propagator 



^,a(w,p) = - 
h 2 



h 2 h 2 A 
8na 47T 2 
0(|q-f|_A)0(| q +E|_A) 9(q-A) 
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At each RG step we first perform the q integration in 
Eq. (17) which is trivial due to the 5- functions in the 
characteristic function \k- Next the frequency integra- 
tion in Eq. (fl7|) is carried out. The computational speed 



is greatly enhanced by mapping the numerical integration 
onto an analytical integration using the spline polynomi- 
als in the interval (— w max , w max ). w max is chosen such 
that the error in the u> integration of the outer regime 
u) > w max introduced due to the approximation ( 16 ) is 



smaller than the accuracy of the numerical solution of the 
system of differential equations (JsJ) . For the final angu- 
lar integration in x = cos 9 on the right-hand side of the 
flow equation (17) we use a numerical integration with 



adaptive nodes in order to cope with discontinuities of 
the integrand. 

In order to obtain a stable numerical result it is suffi- 
cient to calculate the flow of the propagators for roughly 
1500 grid points (u>i,Pj) in frequency and momentum 
space. The corresponding system of ordinary differen- 
tial equations is straightforwardly solved using a Runge- 
Kutta algorithm, which we have implemented in a ver- 
sion with adaptive stepsize in RG time t = ln(fc/A). An 
adaptive stepsize is essential in order to detect the kinks 
in the RG flow due to the sharp Fermi surface of the f- 
fermions at zero temperature. We observe that about 10 4 
RG steps are necessary to obtain an error smaller than 
e ~ 10~ 5 . 

Finally, when the flow reaches the infrared, k = 0, 
we end up with the full Matsubara Green's functions 
G^k=o((jJ,p) and Gr^fc=o(w,p). The initial value of fj,± 
for a given kpa is adjusted such that a vanishing macro- 
scopic occupation of .J,- atoms and molecules is obtained 
at the end of the flow, as discussed in section [TTJ. In order 
to access the spectral functions we perform the analytical 
continuation to real frequencies using a Pade approxima- 
tion. 



We now present our numerical results for the spectral 
functions of the polaron and molecule across the whole 
transition region. First, the Matsubara Green's func- 
tions at the end of the RG flow, Gj.m^oCw, p), are con- 
tinued analytically to retarded Green's functions of real 
frequency, G;/^r(w,p), using the Pade approximation. 
The spectral functions are defined as 



A l/<f{u,v) = 2ImG i/0!R (w,p). 



(19) 



In Fig. [2] the zero-momentum spectral functions 
A|./<*( W !P = 0) are shown as functions of frequency and 
coupling (kpa)^ 1 . 

The coherent single-particle excitations at zero mo- 
mentum are determined by the solutions w qp of the equa- 
tion 



0) 







(20) 



for lu in the complex lower half-plane. Near each quasi- 
particle pole the retarded propagator can be approxi- 
mated by the form 



G;/0, R (w,p = 0) 



J l/<f> 



to 



iO 



(21) 



where the real part of w qp determines the quasi-particle 
energy 



E qp = H + Re[w qp ]. 



(22) 



We have shifted the ground-state energy, which is zero 
in our calculation (vanishing gap), to the conventional 
value Hi from Eq. (12). The imaginary part of the pole 
position determines the decay width 



- Imfcjn 



(23) 



A Fourier transform in time relates the decay width to 
the quasi-particle lifetime 



'qp 



h/r. 



qp- 



(24) 



The quasi-particle weight Zuq is obtained from the fre- 
quency slope at the complex pole position, 



7- 1 



(25) 



Note that an alternative definition of the decay width 
in terms of the self-energy evaluated not at the complex 
pole position but on the real frequency axis, 



all 



Im£|,/^ R (u;,p = 0) 



UJ— Rc UJn 



(26) 



agrees with our definition for T qp only for a single quasi- 
particle pole (21) with Z = 1. However, for the po- 



laron problem there are further excited states and Z < 1. 




FIG. 2: (color online). Spectral functions at zero momentum in dependence on l/(kira). (a), (b): Polaron spectral function 
A±(lj,p — 0). (c), (d): Molecule spectral functions A^(u,p = 0). In order to make the 5-function peak for the ground state 
visible we introduced an artificial width of the quasi-particles of 0.007 cf- 



7 



polarouic side 



molecular side 




■ 5 
0.0 



0.2 0.4 0.6 0.8 

l/{k F a) 



FIG. 3: (color online). Energy spectrum of the single-particle 
excitations of the polaron-to-molecule transition. In the in- 
set we show our fRG result for the ground-state energy (with 
the universal dimer binding energy Eb = —h 2 /(ma 2 ) sub- 
tracted) in comparison to the results obtained with diagMC 
by Prokof'ev and Svistunov [2] (symbols). 



Hence, only r qp from Eq. ( 23 ) can be interpreted as the 



half-width of the peaks in the spectral function and as 
the inverse lifetime. 

We will now in turn discuss the features seen in the 
spectral functions: the peak position (E qp ), width (r„p) 
and weight (Z), first for the polaron (upper row of Fig. [2]) 
and then for the molecule (lower row). 



A. Attractive and repulsive polaron 

Energy spectrum. Let us first look at the energy 
spectrum of the quasi-particle excitations depicted in 
Fig. [3] in dependence on the coupling strength (kpa)^ 1 . 
From our data we find two coherent quasi-particle states 
for the J,-atom, the attractive and the repulsive po- 
laron, and one bound state for the molecule. The at- 
tractive polaron (red solid line) is the ground state for 
(kpa)^ 1 < (fci?a c ) _1 (polaronic side) but becomes an 
excited state for (fcpa)" 1 > (kpac) -1 (molecular side). 
Conversely, the molecule is the ground state on the 
molecular side (blue dashed line) and an excited state 
on the polaronic side, in accordance with the discus- 
sion at the end of section |TlJ For the critical coupling 
strength we obtain (kpa c ) — 0.904(5), which agrees 
with the value (fc^ac)" 1 = 0.90(2) obtained using di- 
agrammatic Monte Carlo (diagMC) by Prokof'ev and 
Svistunov [2]. As shown in the inset of Fig. |3j also the 
values for the energies agree well with diagMC (sym- 
bols). At unitarity, (kpa)~ l = 0, we obtain the ground- 
state energy /xj. = — 0.57ef while diagMC gives the value 
/ij, = — 0.615 ep. Having treated the full frequency and 
momentum dependence of the propagators in the trun- 
cation Q, we can attribute the residual deviation in 
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FIG. 4: (color online). Decay channel for the repulsive po- 
laron. (left) Two-body process which leads to the decay of 
the repulsive polaron. (right) Corresponding contribution to 
the ^-atom self-energy via the optical theorem. 



the ground state energy to the omission of the terms 
^V'^V^V'T and 4>*tp^(j>ip a . The latter term describes the 
atom-dimer interaction and is expected to further reduce 
the ground-state energy in the transition regime in accor- 
dance with the results from the variational wave function 
approach [TO] [TTJ. The term xp^ip^ip^tpf, generated by 
particle-hole fluctuations, is expected to give the main 
correction in the unitarity regime [5]. Both terms can 
be included in the fRG flow as additional flowing cou- 
plings, or implicitly by using the Katanin scheme |40| or 
rebosonization |41) . 

Until recently [53] most experiments with ultracold 
Fermi gases have focused on the lower, attractive branch 
on the BEC side (kpa)^ 1 > 0. There exists, however, 
also the repulsive polaron branch (solid green line) which 
corresponds to a higher excited state of the ^-atom inter- 
acting repulsively with the f-Fermi sea. Our results for 
the energy of the repulsive branch agree with the weak- 
coupling results [IS] for (fcpa) -1 > 1. In the strong- 
coupling regime our energies lie between the result from 
the non-self-consistent T-matrix approach [21j and the 
MC results for square well potentials 24 . In the polaron 
spectral function, cf. Fig. [2]^h), one can clearly discern 
the attractive polaron branch as a very sharp peak at 
low frequencies, and the much broader repulsive polaron 
branch at higher frequencies. 

Decay widths. The repulsive polaron has a large 
decay width r rep , as calculated from Eq. (23) and de- 
picted in Fig. [5] and correspondingly 



a short lifetime. 
The leading-order decay channel for the repulsive po- 
laron is the process shown in Fig. Qleft) where the re- 
pulsive polaron, which is an excited state, decays to the 
attractive, and energetically lower lying, polaron due to 
the interaction with an f-atom. This diagram can be 
translated via the optical theorem into a contribution 
to the imaginary part of the |-atom self-energy, as de- 
picted in Fig. fright). This self-energy diagram is al- 
ready included in the non-self-consistent T-matrix prop- 
agator P£ sc , and has been studied recently using this ap- 
proximation [5T] . Of course, this diagram is also included 
in our fRG approach. 

In the weak-coupling limit (kpa)~ l — > 00 the excita- 
tion becomes sharp, r rcp — > 0, and the repulsive polaron 
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FIG. 5: Decay width r rep of the repulsive polaron as a func- 
tion of the coupling (^a) -1 - We also show the width accord- 
ing to the approximate formula Eq. (261 (dotted line). Inset: 
Decay width r a tt of the attractive polaron. 



is a well-defined quasi-particle. Toward unitarity, r rcp 
grows but remains a well-defined, finite quantity even at 
unitarity. We find that indeed r rop , and not the approx- 
imation T^„ from Eq. (26), is the correct half-width at 



rep 



half-height of the respective peak in the polaron spectral 
function in Fig.^lb). For (kpa)^ 1 < 0.6 the energy E rep 
of the repulsive Dranch exceeds the bath Fermi energy, 
-E'rep > £f- At this point it is energetically favorable to 
spin-flip the impurity atom, which can be interpreted as 
the condition for the onset of saturated ferromagnetism 
[2"0ll26| . At the same time the decay width r rop > 0.2 ep 
is large, which potentially destabilizes a ferromagnetic 
phase [25] . 

On the polaronic side (kp a)" 1 < (ifOc)" 1 the attrac- 
tive polaron is the stable ground state with decay width 
T a tt — 0, while on the molecular side it is an excited 
state with finite lifetime and decay width r at t > 0, see 
Fig. [^inset). This decay is much weaker and also quali- 
tatively different from the repulsive channel. The attrac- 
tive polaron can decay by a three-body recombination 
process as shown in Fig. [6](left) . Via the optical theo- 
rem this process can be translated into a contribution to 
the |-atom self-energy as depicted in Fig. fright), plus 
an additional contribution with crossed lines. This de- 
cay channel has recently been studied using an explicit 
three-loop calculation [17]. The resulting finite lifetime 
cannot be seen in the non-self-consistent T-matrix ap- 
proximation, where the self-energy corrections of the J.- 
atom — indicated by the inner white box — are not fed 
back into the T-matrix ~ P^ 1 [5] [5Tj. In contrast to 
the non-self-consistent T-matrix calculation, our fRG in- 
cludes the full feedback of both the I and <f> self-energies, 
denoted by bold internal lines in the flow equations in 
Fig. [T] Therefore, the contributions from the decay dia- 
gram in Fig. fright), and many more, are automatically 
included in our approach. 

Quasi-particle weights. Fig. [7] depicts the quasi- 
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FIG. 6: (color online). Decay channel for the attractive po- 
laron. (left) Three-body recombination process which leads 
to the decay of the attractive polaron. (right) Corresponding 
contribution to the 4,-atom self-energy via the optical theorem 
(there is also a contribution with crossed lines). 



particle weights of the attractive and repulsive polaron 
computed using Eq. ( 25 ) . On the polaronic side the at- 



tractive polaron state contains most of the weight, but 
as one moves toward the molecular side the spectral 
weight gradually shifts to the repulsive branch, and the 
corresponding peak in the polaron spectral function in 
Fig. [2jb) becomes larger. We find that the attractive and 
repulsive branches almost completely make up the total 
spectral weight, hence the contribution from the incoher- 
ent background is very small. This is also apparent in the 
polaron spectral function in Fig. [2^b). 

Our results for the quasi-particle weights agree well 
with those from the non-self-consistent T-matrix and 
variational wave-function approaches. At unitarity we 
obtain for the attractive polaron 2j. ia tt = 0.796 compared 
to the variational value Z^att — 0.78 OHO]- For the re- 
pulsive polaron at (kpa)^ 1 = 1 we find Zi >Iep — 0.71, in 
agreement with the recent non-self-consistcnt T-matrix 
calculation [21] . 

Note that there is an alternative definition of the quasi- 
particle weight [TU] 



= lim G^,p = 0) . 



(27) 



This definition has to be treated with care: on the molec- 
ular side of the transition the polaron acquires a finite 
decay width > but nonetheless continues to be a 
well-defined quasi-particle with finite spectral weight, as 
can be seen from Fig. [2} However, definition (27) yields 



zero as soon as > and in this case cannot be in- 
terpreted as a measure of spectral weight anymore. In 
contrast, the definition (25) remains correct for a finite 



decay width and accordingly our data for Z± shows no 
discontinuity at the transition. In the experiment the fi- 
nite lifetime of the attractive polaron on the molecular 
side complicates the direct measurement of by radio- 
frequency spectroscopy because the molecular state and 
not the attractive polaron becomes occupied as the initial 
state (cf. section [v|. 
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FIG. 7: (color online). Quasi-particle weight of the at- 
tractive (solid line) and repulsive (dashed line) polaron. The 
weights of the two quasi-particle peaks in the J, spectral func- 
tion almost completely make up the total spectral weight, 
and the contribution from the incoherent background is very 
small. 



B. Molecule 

Energy spectrum. The molecule spectral function 
in Fig. [2^d) displays a sharp quasi-particle peak of the 
bound state at low frequencies, followed by an incoherent 
background at higher frequencies which actually carries 
most of the spectral weight. Note that this background is 
not taken into account in the simple derivative expansion 
in the appendix nor in the Wilsonian RG approach [T5] , 
On the molecular side (fc^a) -1 > (kp a c ) _1 the molecule 
is the ground state and is clearly separated from the inco- 
herent continuum. On the polaronic side of the transition 
the molecule becomes an unstable, excited state and de- 
velops a clearly visible finite decay width in the spectral 
function. 

Decay widths. The leading decay channel of the ex- 
cited molecule state is via the three-body recombination 
process shown in Fig. [8](lcft) . Via the optical theorem 
this process can be translated into a contribution to the 
molecule self-energy as depicted in Fig. fright). Simi- 
larly to the attractive polaron, in the non-self-consistent 
T-matrix approximation the J.-atom self-energy correc- 
tions in are not fed back into the T-matrix ~ P, 



and the molecule does not decay. In contrast, the dia- 
gram Fig. fright) is included in the fRG, which leads to 
the visible broadening in the spectral function. 

Ref. [T7j shows by an analytical calculation of the phase 
space for three-loop diagrams of the type in Fig. [8] that 
the decay width of the molecule scales as 



r„ cx Aw 9 / 2 



Aw = Ea 



E 



(28) 



where Aw is the difference between the energy levels of 
the excited molecule and the attractive polaron ground 
state. In Fig. [9] we show as a function of Aw in a dou- 
ble logarithmic plot. The large fluctuations of our nu- 
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FIG. 8: (color online). Leading decay channel for the ex- 
cited molecular state, (left) Three-body recombination pro- 
cess which leads to the decay of the molecular state, (right) 
Corresponding contribution to the molecule self-energy via 
the optical theorem (there is also a contribution with crossed 
lines). 



merical data are due to the accuracy of the Runge-Kutta 
integration as well as due to the restriction to a finite 
number of Matsubara frequencies. We have estimated 
the corresponding error by comparing the results for dif- 
ferent grids with varying number and position of (Mat- 
subara) frequencies. The solid line in Figj9] indicates the 
power law Aw 9//2 . The triangles in Fig. |9jcorrespond to 
a calculation with a higher number of frequencies and we 
find a convergence to the solid curve for decay widths 
larger than our numerical integration accuracy e = 1CP 5 . 
This indicates that the error for larger can be at- 
tributed to the Pade approximation, while for < e the 
accuracy of our results becomes limited by the absolute 
error of our numerical integration. With our fRG calcu- 
lation we are thus able to verify the prediction by Bruun 
and Massignan [17] , and the correctness of the power law 
attests to the strength of our method to describe many 
features of the polaron-to-molecule transition in one uni- 
fied approach. 

Quasi-particle weight and compositeness. The 
quasi-particle weight Z^ of the molecular bound state 
in the spectral function in Fig. u\d) is very small, Z^ w 
0.002 at unitarity, and increases slowly toward the molec- 
ular limit, see Fig. 10 Inset). For broad resonances, 
h 2 ~ AB — > oo, the two-channel model ^ is equiva- 



lent to the single-channel model ([I]) where Z^ — [3T]. 
Specifically, we obtain for the weight of the bound state 
in vacuum 



Z<h = 



32vr 



(vacuum, A — > oo). 



(29) 



This is consistent with the interpretation of Z§ as the 
closed-channel admixture (cf. Eq. (29) in [S]). In our 
calculation we set the physical UV cutoff scale to A = 
10 3 k F , which is of the order of the inverse Bohr radius, 
and choose h 2 < oo. We observe Z$ oc l/(kpa) on the 
BEC side and a deviation from the vacuum scaling close 
to unitarity, which may be due to a combination of finite 
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FIG. 9: The decay width F^> of the excited molecular state as 
a function of the energy difference Alo — E$ — -E^att between 
the excited molecule and attractive polaron ground state. The 
solid line indicates the power-law scaling oc Aw 9 ' 2 . 
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FIG. 10: Compositeness (1—Z$) of the molecular bound state. 
A value of 100% would indicate that the molecule has no over- 
lap with elementary closed-channel bosons. Inset: Molecular 
residue 



V. RF-RESPONSE OF THE 6 LI FERMI GAS 



density corrections and the admixture of closed-channel 
molecules in the microscopic action by choosing finite 
values of h and A. 

Furthermore, the quasi-particle weight Z can be inter- 
preted as the overlap between the "true" particles and 
the "elementary" , or bare, particles in the microscopic 
action The attractive Fermi polaron becomes ele- 
mentary, Z^. a tt - > 1, in the BCS limit {kpa — > _ ), while 
in the opposite limit of kpa — > + the repulsive Fermi 
polaron becomes elementary, 2j.,rep — > 1- Near unitarity, 
both excitations have a sizable weight. In contrast, the 
molecular bound state is almost exclusively a composite 
particle in the whole transition region. Indeed, the de- 
viation of the quasi-particle weight from unity, 1 — Z^, 
is a well-established measure of compositeness in nuclear 
physics [3D] , and in Fig. 10 we show that for our choice of 



the Yukawa coupling h the compositeness of the molecule 
is very large (> 98%). This is consistent with the mea- 
surement of a small molecular weight for a balanced 
6 Li Fermi gas close to a broad Feshbach resonance by Par- 
tridge et al. |42j . In experiments with a narrow Feshbach 
resonance the compositeness will decrease and a single- 
channel description becomes invalid. A strength of the 
fRG approach is that both situations are naturally de- 
scribed by tuning the values of h and G^a. 

In the vacuum there exists no molecular state for neg- 
ative scattering length a as can easily seen from Eq. Q , 
which has no bound-state pole. In the presence of a 
medium of f-fermions, however, the molecule propagator 
develops an excited bound-state pole also for negative 
scattering length a. This bound state will become the 
superfluid ground state if the impurity density exceeds a 
critical threshold [TU] , 



As a final application we want to connect our results for 
the polaronic spectral function to experimentally observ- 
able radio-frequency (rf) spectra. The attractive branch 
of the polaron-to-molecule transition has been studied 
experimentally by Schirotzek et al. [6] using a popula- 
tion imbalanced, two-component mixture of 6 Li atoms. 
In the experiment the rf-response of the system has been 
used to infer information about the low-frequency behav- 
ior of the fcrmionic spectral functions. For instance, the 
ground-state energy and the residue Z± of the 4,-fermions 
were measured and confirmed the theoretical predictions. 

In order to measure the rf-response, an rf-pulse is ap- 
plied to the system which drives the transition of the 



atoms to a third, initially empty state. In Fig. 11 we show 
the scattering length profile of 6 Li versus the magnetic 
field. In the experiment [B] a mixture of fermions initially 
in the hyperfine states |1) and |3) had been prepared in 
a range of the external magnetic field B = 630 . . . 690 G 
(shaded area). In this regime the scattering length a 13 
in the initial state is large and positive, while the final 
state scattering lengths an, ^23 are rather small. 

Information about the spectral function can be ac- 
cessed from rf-spectroscopy, for example by populating 
the particle under investigation up to the energies one is 
interested in and then driving the transition to a weakly 
interacting final state for which the spectral function is 
well-known. This route had been taken for the study 
of the attractive polaron. The repulsive polaron branch, 
on the other hand, has not yet been observed directly 
in experiments. The main complication is that the re- 
pulsive polaron has a very short lifetime in the strong- 
coupling regime of interest (cf. Fig. [5]). Hence, its macro- 
scopic population is inhibited on longer time scales, and 
even after a fast ramp to the desired magnetic field most 
minority atoms will have decayed into the respective 
ground state. A similar situation arises in the detection 
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FIG. 11: (color online). Scattering length profile of 6 Li atoms 
in the three lowest hyperfine states as calculated by P. S. 
Julienne based on the model described in [43]. The stars 
indicate the magnetic field values for which we determine the 
rf-response of the imbalanced Fermi gas in Fig. |13| 



of Efimov trimers in a three-component mixture of 6 Li 
atoms [HI US] . The decay of the repulsive branch is also 
of relevance for the balanced system and the competi- 
tion between ferromagnetic order and molecule formation 

In this section we propose an experimental procedure 
to circumvent these difficulties and directly analyze the 
spectral function of the repulsive polaron. A strongly 
imbalanced two-component 6 Li Fermi gas is prepared in 
hyperfine states |1) and |3) for magnetic fields B > 690 G 
across the (1,3) Feshbach resonance. In this regime the 
initial scattering length is negative, 013 < 0. One then 
drives an rf-transition to the final state |2) which is char- 
acterized by large, positive scattering lengths 012 and 
(Z23 and thus strong interactions. Because the attrac- 
tive polaron spectral function of the initial state, with its 
negative scattering length 013, is well understood both 
experimentally and theoretically, the final-state spectral 
function can then be analyzed in a controlled fashion. 

Within linear response theory the induced transition 
rate from the initial state |i) to the final state |/) is given 
by 



I(u) L ) = 2fZ 2 Im xr(m/ - M» ~ <*>i) 



(30) 



where the Rabi frequency f2 is given by the coupling 
strength of the rf-photon to the atomic transition, fii 
(pf) the initial (final) state chemical potential, and wl 
denotes the rf-frequency offset with respect to the free rf- 
transition frequency. Neglecting the momentum of the rf- 
photon, the retarded rf-susceptibility xr is given by the 
analytical continuation to real frequencies of the correla- 
tion function in imaginary time r (Matsubara frequency 



e^(T T ^(r,r)^(r,r) 

x^(r',0)V/(r',0)), (31) 





B field [G] 


(fc F ai 3 ) _1 


(k F a,23)~ 1 


(fc F ai 2 ) _1 


(a) 


810.3 


-1.88 


0.0 


0.39 


(b) 


800.8 


-1.80 


0.2 


0.58 


(c) 


788.2 


-1.70 


0.5 


0.86 


(d) 


780.6 


-1.62 


0.7 


1.04 



TABLE II: Int erac tion parameters at the four transitions in- 
dicated in Fig. 11 using kpf = 0.00015 Oq~ . 



where 
Eq, 



T T is the imaginary time-ordering operator. 
(31) leads to various diagrammatic contributions 



which are in general difficult to handle if the final-state 
interactions are not negligible [49]. Here we will calcu- 
late Eq. (31 1 in a simple approximation with full Green's 



functions but without vertex corrections. In this approx- 
imation Eq. (31) yields the susceptibility in Matsubara 
frequency 

X(«) = / Gi(k,i/)G/(k,i/ + w). (32) 
The rf-response in real frequency is then given by 



/K) = n 2 



Ik JO 27T 

x A f (k,v) Ai(k,u + fi f 



Hi - uj l ) (33) 



where \ij — jj,f + ujl > 0- 

In Eq. (33) the initial-state spectral function A; is 
probed for negative frequencies only. Hence, there is 
no rf-response for the pure polaron problem at vanish- 
ing density and chemical potential /ij '*- In the experi- 
ment one has, however, a small but finite concentration 
x = n^/n^ of J,-fermions which leads to an observable 
rf-response. In order to describe the experimental sit- 
uation we therefore need a calculation for a finite mi- 
nority (1) density characterized by a chemical potential 
> (i^ . Within the fRG framework such a 



calculation requires the regulator R± in Eq. ( 10 1 to be 
adjusted in order to cope with the finite Fermi surface 
of ^-fcrmions which complicates the computation. For- 
tunately, our calculation for the polaron problem shows 
that for negative scattering length, where the polaron 
is the ground state and decay processes do not matter, 
the fRG results are in excellent agreement with the re- 
sults from a non-self-consistent T-matrix approach. We 
may therefore use this approach instead of a full-feedback 
fRG for the calculation of the imbalanced Fermi gas of 
finite densities and n_\ in order to determine the initial- 



state spectral function. The non-self-consistent calcula- 
tion was also done by Punk and Zwerger [46] and is ob- 
tained in our fRG formulation by simply switching off 
the feedback of the 4,-atom self-energy into the molecule 
flow. 

Because the occupation of 4,-atoms is small only the 
low- momentum modes are relevant. The 4,-atoms form a 
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FIG. 12: (color online) . Full momentum and frequency depen- 
dence of the polaron spectral function Ai(uj,p) at unitarity 



degenerate Fermi gas of polaronic quasi-particles and the 
spectral function can be approximated by [TB] 



Ai(ui, p) = 27rZi (51a; — 



2m 



A 



-I 



(34) 



where Zi is the residue and m| is the effective mass of 
the impurity atoms. A determines the impurity concen- 
tration x. 

We have calculated the parameters ml, Z i} and A 
as functions of (kpa)^ 1 via the non-feedback (non-self- 
consistent T-matrix) calculation, fii is determined self- 
consistently to ensure the correct impurity density. In- 



serting the spectral function Eq. (34 1 into the suscepti- 



bility ( 33 ) we obtain for the rf- response 



n 2 z, 



V 2m l A 



2ir 2 



Mi 



dkk 2 



H f + lo l - A 



2mt 



(35) 



For the final state |/) we use the spectral function ob- 
tained in section IIVI from the full fRG calculation. In 
Fig. [12] we show an example of the full final-state spec- 
tral function in dependence of frequency and momentum 
which enters the momentum sum in Eq. (35). One can 



clearly discern the broadening of the attractive polaron 
branch at larger momenta, as well as the broad repulsive 
branch at higher frequencies. A similar feature appears in 
the spectral function for the balanced Fermi gas above T c 
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FIG. 13: (color online), rf-spectra I{ijJl) for minority species 
|1) (solid blue lines) and minority species |3) (dashed red 
lines). The interaction parameters correspond to the four 
magnetic field values marked in Fig. |ll| and listed in Table \U\ 
(a) top — (d) bottom. 



as calculated by Haussmann et al. l3H]. The fermionic 
spectral functions and rf-spectra in this work were de- 
termined using state-of-the-art self- consistent T-matrix 
(2PI) approximations and form the basis for a full linear 
response calculation including vertex corrections [5U] , 

Initially, the gas is prepared in a |1) and |3) mixture, 
and both states can serve as minority or majority species, 
for example, |f) = |3) and \\) = |1) such that \i) = |1) 
and I/) = |2). The initial occupation of the minority 
states is small and for our numerical calculation we use 
x = ni/rif = 0.01. The energy scale is set by the Fermi 
momentum of the majority species kpf — 0.00015 a^ 1 as 
appropriate for the MIT experiment [5] . We calculate the 
rf-spectra at four values of the magnetic field indicated 
as stars in Fig. |TTJ and listed in Table [TTJ. The resulting 
spectra are shown in Fig. |13| the solid blue lines indicate 
the response for minority species |1), while the dashed red 
lines correspond to minority species |3). In the latter case 
the sign of the frequency offset ujl is changed because the 
state 1 3) is energetically above the final state. 

The position of the sharp attractive polaron peak at 
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negative frequency offset ljl shifts with the interaction 
parameter (kpaf)^ 1 in accordance with the energy spec- 
trum Fig. [3] One observes that the attractive polaron 
loses quasi-particle weight on the molecular side, cf. 
Fig. [7J In contrast, the repulsive polaron branch gains 
quasi-particle weight toward the molecular side, and the 
respective peak in the rf-spectra becomes both larger and 
narrower, and one can read of the increasing lifetime. 

The spectra in Fig. [13] are convolved with a sine func- 
tion sinc 2 (wT/2) which gives the response to an rf-pulse 
with a rectangular profile of length T = 20 ms [35] . While 
our curves are computed for zero temperature, a finite 
temperature ~ 0.01 T F would lead only to a slight broad- 
ening of the experimental rf-peaks. The broadening of 
the attractive polaron due to the finite lifetime on the 
molecular side (kpaf)^ 1 > (fci?a c ) _ , however, is too 
small to be resolved. Note that we have included both 
final-state and initial-state interaction in our calculation: 
the knowledge of the spectral function for the initial state 
allows for a detailed study of the final-state spectral func- 
tion. 



VI. DISCUSSION 

We have presented a new computational method to 
solve the non-perturbative, exact renormalization group 
equation ([3]) and have demonstrated its efficiency for the 
Fermi polaron problem as a specific example. The inclu- 
sion of the full frequency and momentum dependence of 
the propagators opens up new perspectives to apply the 
functional renormalization group to problems where the 
detailed dynamics of the relevant degrees of freedom be- 
comes important [51j . In particular, the method draws 
its strength from the possibility to successively bosonize 
further channels of the interaction via additional aux- 
iliary fields (Hubbard-Stratonovich transformation) [52] . 
In this way one can partially capture the complicated an- 
alytical structure of higher-order vertex functions r'"', 
including possible quasi-particle poles and branch cuts, 
as we have explicitely shown for the s-wave scattering 
channel in the polaron problem. In combination with 
the recently developed flowing rebosonization technique 
[41] our numerical method can be extended to also incor- 
porate re-emerging vertices. Our approach complements 
the proposal for bosons [29] and additionally includes 
fermions. 

For the Fermi polaron problem we achieve a unified 
description of many dynamical effects beyond thermo- 
dynamics. We verify the non-trivial power-law scaling 
of decay rates [T7] and determine the properties of the 
repulsive polaron with a method beyond the non-self- 
consistent T-matrix approximation |201 121) . This is of 
value in the ongoing debate about the possible occur- 
rence of ferromagnetism in ultracold Fermi gases with 
short-range interactions [20T426j . The polaron problem 
sheds light on this question in the limit of strong popu- 
lation imbalance. For the repulsive polaron we find the 



critical interactions strength kpa = 1.57 from our numer- 
ical data (Fig. [3]) . Going to a finite density of 4,-fermions 
is straightforward within the fRG and involves only a 
slight modification of the regulator of the ^-fermions ( 10 ) 



as long as no spontaneous symmetry breaking occurs. By 
continuity we can infer that the repulsive branch will re- 
main to exist for small but finite .^-population and will ex- 
ceed the critical energy cf for the presumed onset of satu- 
rated ferromagnetism. It is an open question whether for 
larger impurity concentrations the repulsive branch is so 
strongly renormalized that saturated ferromagnetism can 
be ruled out |26) . or whether competition with molecule 
formation may preclude the observation of ferromagnetic 
domains |25j . Answering these questions will require a 
full non-equilibrium calculation. 

There has been much theoretical progress on the re- 
pulsive Fermi gas with short-range interactions, but rel- 
atively few experiments have been completed. While the 
repulsive 3 He Fermi gas has been studied extensively in 
experiment, it is not dilute and has a large repulsive hard- 
core potential [53 . In contrast, ultracold Fermi gases 
offer the realization of a proper contact interaction of 
tunable strength. We predict rf-transition rates for the 
repulsive branch and propose a possible route to measure 
these excited states in a 6 Li Fermi gas. This is a chal- 
lenging problem because the final repulsive polaron state 
is highly unstable. One possible approach could involve 
fast tomographic imaging similar to the MIT experiment 
[6] , another might be to measure the loss in the final state 
which is expected to scale with the rf-transition rate for 
a constant rf-pulse time but may be suppressed by the 
quantum Zeno effect [U]. Hence, the possible observa- 
tion of the repulsive polaron represents not only a test of 
theoretical predictions but poses an interesting challenge 
touching several aspects of many-body physics. 
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Appendix A: Derivative expansion 

In this appendix we qualitatively study the polaron- 
to-molecule transition using a simple approximation for 
the renormalized Green's functions Gj, £ and G^^. As we 
have seen, G^k and G^k generally develop a complicated 
frequency and momentum (w,p) dependence. Here, we 
will use an expansion in small frequencies and momenta 
(derivative or gradient expansion) which allows for an an- 
alytical evaluation of the loop integrals on the right-hand 
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side of the flow equation (JSj) _ The derivative expansion 
proves to be a good approximation if one is interested in 
physics determined by the structure of the Green's func- 
tions close to their poles as, for example, in the descrip- 
tion of phase transitions and critical phenomena |32| . In 
the derivative expansion the dependence of the inverse 
propagators Pjm^u;, p) on the RG scale k is approxi- 
mated by 



Pl,k(u, p) = A i>k • (-iw + p 2 ) + m 2 k 
P^fc(w, p) = A^ k • (-iw + p 2 /2) + ml k . 



(Al) 



propagator Q. We therefore focus only on the correct 
calculation of the scattering length a for q = 0, which 
leads to the infrared condition («i^ a jr, =0 ) 2 = —h 2 /(8ira) 
for /ij, = /i-j- = 0. In this case of two-body physics P\.,k 
is not renormalized, A^k = 1 and wi 2 fc = 0, and the dif- 
ferential equations for A^k an d ni^ k decouple and can 



be solved analytically. The integration of ( A2 1 in the 
vacuum limit yields 



(in 



vac \2 
<p,k=0J 



h 2 A 
4^' 



(A3) 



In this approximation, which is similar to the one used 
in the Wilsonian RG approach to the polaron problem 
[15j , we assume that the renormalization of the frequency 
and momentum coefficients is given by common wave- 
function renormalizations A^k and A^, t k, respectively, 
which in turn are related to the quasi-particle weights 



via Z l/(f> = A 



i/<t>' 



The flowing gap terms m| k and 



4>,k 



arc related to the flowing static self-energy via Eq. ( 11 ). 

From equation (JsJ) one can derive the flow equations 
for the four running couplings A^ k , A^k-, m lk> an d 



l t,fc 



The f-fermions are not renormalized, A^ t k — 1 an d 



-//-)-. With the sharp cutoff (10), the frequency 



as well as the momentum integrations can be performed 
analytically. The resulting flow equations read 



d k Ax = - 



2h 2 k 



6(ny - 2k 2 ) 



- fc2 



(k 2 + Mt 



(-fc 2 + 2/i t + 2m 2 /A^f 



d k m 2 



h 2 k 
ir 2 A d 



2k 2 ) 
k 



V 2m 2 M ) 2 

(A2) 

It 2 



k 2 + ^ + 2m 2 /At 



-k 2 



dkA^ = - 



dkm 2 , 



h 2 k 
2tt 2 A 

h 2 k 



2fPf + 2m 2 /A^ 



k 2 

(2k 2 +^+m 2 /A i f 



2w 2 A l (2k 2 + ^ + m 2 /A l )- 



The initial conditions for this system of differential equa- 
tions (A2) are specified at the UV scale k = A. We 



note that the derivative expansion (Al) of the molecule 



propagator P^^(a;, p) cannot account for the correct vac- 
uum scattering amplitude (p)]), because the term iq is 
rooted in the non-analytical structure of the molecule 



This leads to the UV condition for the molecule gap, 



8tt 



(2A/tt - a" 1 ), 



(A4) 



which incorporates the correct regularization of the UV 
divergence A in Eq. (A3). At the UV scale the momen- 



tum and frequency dependence of P^,,a can be neglected 
due to the large bosonic gap m 2 , A , and we set A 



1. 



For a finite de nsit y of f-atoms the system of differ- 
ential equations ( A2 ) is solved numerically. The initial 



values for the 1 propagator are then given by m 2 A = — /ij, 
and Aua = 1. The down chemical potential /xj, is deter- 
mined in the way discussed in section [TT] We find that 
the polaron is indeed the ground state for interaction 
strengths (kpa)^ 1 < (kpac)^ 1 whereas the molecule be- 
comes the ground state for (kpa)~ l > (kpa c ) , with 
{kpa,)- 1 = 0.96. 

The energy spectrum from the simple derivative ex- 
pansion is in qualitative agreement and even in rough 
quantitative agreement with the results obtained from 
our new numerical method and other theoretical calcula- 
tions. A drawback of the derivative expansion is that 
it is impossible to extract a reasonable spectral func- 



tion from an ansatz of the form ( Al ) as it only accounts 
for a single coherent quasi-particle excitation. Neither 
higher excited states, such as the repulsive polaron, nor 
the incoherent background, which comprises the major 
weight for the molecule (cf. section IV), can be cap- 
tured with this ansatz. Furthermore, although decay 
processes lead to a finite lifetime of the excited polaron 



and molecule branches, in the simple approximation ( Al ) 



these states have a vanishing decay width. Note that the 
finite lifetime of excited states can be obtained neither 
from simple variational wave functions [7] nor from the 
non-self-consistent T- matrix approach [8] . Their descrip- 
tion requires the full self-energy feedback developed in 
section Mil 
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